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ABSTRACT 

We present optical HST/STIS spectroscopy of RZ 2109, a globular cluster in the elliptical galaxy 
NGC 4472. This globular cluster is notable for hosting an ultraluminous X-ray source as well as 
associated strong and broad [O ill] AA4959, 5007 emission. We show that the HST/STIS spectroscopy 
spatially resolves the [O ill] emission in RZ 2109. While we are unable to make a precise determination 
of the morphology of the emission line nebula, the best fitting models all require that the [O III] A5007 
emission has a half light radius in the range 3-7 pc. The extended nature of the [O III] A5007 emission 
is inconsistent with published models that invoke an intermediate mass black hole origin. It is also 
inconsistent with the ionization of ejecta from a nova in the cluster. The spatial scale of the nebula 
could be produced via the photoionization of a strong wind driven from a stellar mass black hole 
accreting at roughly its Eddington rate. 

Subject headings: galaxies: individual: NGC 4472 - galaxies: star clusters - globular clusters: general 
- star clusters: individual: RZ 2109 - X-rays: binaries - X-rays: galaxies: clusters 



1. INTRODUCTION 

Some of the best evidence for black holes (BHs) in glob- 
ular clusters (GCs) has come from X-ray studies of extra- 
galactic GCs. Specifically, X-ray sources with luminosi- 
ties in excess of the Eddington limit for an accreting neu- 
tron star have been observed - suggesting accretion on to 
a more massive BH primary (e.g., Sarazin et al. 2000; An- 
gelini et al. 2001; Fabbiano 2006). Some of these ultralu- 
minous X-ray sources are observed to vary on short time 
scales, which is important because it eliminates the pos- 
sibility that the high X-ray luminosity arises from the su- 
perposition of multiple accreting neutron stars (Kalogera 
et al. 2004). To date, six GC X-ray sources have been 
proposed as accreting BHs based on these two key prop- 
erties (luminosity and variability). These sources reside 
in GCs in the galaxies NGC 4472 (two systems; Mac- 
carone et al. 2007, 2011), NGC 1399 (Shih et al. 2010), 
NGC 3379 (two systems; Brassington et al. 2010, 2012) 
and NGC 4278 (Brassington et al. 2012). In addition 
to these objects, BHs have been proposed to exist in 
GCs in two other galaxies. A BH candidate was identi- 
fied in a GC in the elliptical galaxy NGC 1399 based on 
non-variable ultraluminous X-ray emission and nebular 
emission lines (Irwin et al. 2010). Nine GCs in M31 have 
fainter X-ray emission, but spectra that look like BH low 
mass X-ray binaries (LMXBs), rather than neutron star 
LMXBs (Barnard et al. 2008, 2012). 

^Based on observations made with the NASA/ESA Hubble 
Space Telescope, obtained from the data archive at the Space 
Telescope Science Institute. STScI is operated by the Associa- 
tion of Universities for Research in Astronomy, Inc. under NASA 
contract NAS 5-26555 



The ultraluminous and variable X-ray source 
XMMU122939. 7-^075333, located in the GC RZ 2109 
in the elliptical galaxy NGC 4472, is a very strong 
BH candidate (Maccarone et al. 2007). In a 2004 
XMM Newton observation, the peak luminosity of 
this source was Lx — 4x lO'^^erg/s, which is an order 
of magnitude greater than the Eddington limit for 
an accreting neutron star. Furthermore, during this 
observation the luminosity varied by a factor 7 in only a 
few hours (Maccarone et al. 2007). The optical spectrum 
of this cluster revealed extraordinary [O III] emission 
(Zcpf et al. 2007). Follow-up observations confirmed 
the presence of extremely broad [O III] AA4959, 5007 
emission lines with velocity widths of ^1500 kms~^ and 
a luminosity of L([0 III] A5007) = 1.4x10^^ ergs/s (Zepf 
et al. 2008). In addition to the extreme luminosity and 
velocity width of this emission, the non-detection of H/3 
from the cluster requires that the nebula is extremely 
hydrogen deficient, with L([0 iii])/L(H/3) > 30. Strong 
forbidden lines are extremely rare in GCs (e.g. Peacock 
et al. 2012) and the emission is very likely related to 
the (equally rare) ultraluminous X-ray source in the 
cluster. As noted by Zepf et al. (2007), this nebula can 
therefore provide important clues as to the nature of 
the accreting compact object in RZ 2109. In this paper, 
we present spatially resolved Hubble Space Telescope 
(HST)/ Space Telescope Imaging Spectrograph (STIS) 
longslit spectroscopy of RZ 2109 (Section 2). In Section 
4, we use these data to investigate the size of the 
emitting nebula and, in Section 5, we compare this with 
the expectations from the current models that have been 
proposed to explain the X-ray and optical emission from 
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RZ 2109. 

2. STIS SPECTROSCOPY 

In order to study the spatial extent of the [O ill] emis- 
sion from RZ 2109, we obtained HST STIS longslit spec- 
troscopy (proposal ID 11703, PI Zepf). The cluster 
was observed on four dates: 2009/12/29; 2010/01/01; 
2010/03/03; 2010/04/11. AU observations were obtained 
through the 52x0.2" slit using the G430L grating (with 
a central wavelength of 4300A). On each date 4 obser- 
vations were taken, resulting in 16 exposures and a total 
exposure time of 41.2 ksec. 

We used the calibrated 2D rectified spectral images 
produced by the standard STIS pipeline routines run 
by CALSTIS (these are the images given 'sx2' suffix by 
the pipeline). This performs the bias subtraction, cos- 
mic ray rejection, dark subtraction and applies the flat 
fields. It then applies a geometric, wavelength and flux 
calibration. The resulting image has a scale of 1 pixel — 
0.05078" in the spatial direction and is wavelength cali- 
brated in the dispersion direction, with 1 pixel = 2. 73 A 
and an estimated line spread function of 1.5 pixels (4.lA). 

The final reduced spectral image was found to have 
significant numbers of hot pixels remaining. We (in com- 
bination with the STIS helpdesk) experimented with re- 
running the CALSTIS pipeline using darks taken with sim- 
ilar temperatures and within the same annealing period, 
but found no improvements over using the master darks 
for our observations. While the darks vary significantly 
with time, the greater S/N of the master dark gave bet- 
ter results. We were therefore unable to correct for these 
pixels. Instead, we flagged the bad pixels by running a 
5cr rejection script over 10 pixel regions in the dispersion 
direction. 

The 16 resulting 2D spectral images were combined 
using the iraf task imcombine. The spatial and wave- 
length calibrations were found to be accurate enough to 
combine the reduced images without applying any shifts. 
The images were scaled and weighted based on their ex- 
posure times and combined together using mask images 
to reject bad pixels. The following analysis is based on 
this final combined image. 

The images are flux calibrated by the pipeline 
to give the surface brightness per A in units of 
erg/cm^/s/A/arcsec^. To produce the ID spectrum of 
RZ 2109, we integrated the spectrum over a region 7 
pixels wide (as is recommended for STIS observations). 
The resulting fluxes were multiplied by the DIFF2PT 
keyword (produced by the pipeline) to give the flux cal- 
ibrated spectrum, corrected for slit losses (assuming a 
point source). Given that our source is slightly ex- 
tended, we experimented with wider extraction regions, 
but found the default resulted in the best S /N spectrum. 
Finally, the fluxes were corrected for known charge trans- 
fer efficiency (GTE) effects (see e.g. Goudfrooij et al. 
2006). This signal loss is significant for these observa- 
tions of a faint target, because they were taken before 
the extent of the GTE problem was commonly known 
and hence centered on the middle of the detector (rather 
than the now commonly used 'El' aperture, located close 
to the detector readout side). The GTE correction was 
calculated as 0.52, based on the equation defined in Sec- 
tion 3.4.6 of the "STIS data handbook" (v6.0, Bostroem 
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Fig. 1.— The STIS spectrum of RZ 2109 (red). The blue line 
shows a previously obtained Keck spectrum of this cluster. The 
black line is a model of a 12 Gyr star cluster with the magnitude 
and metallicity of RZ 2109. The zoomed panel shows the 5007A 
emission line region (located at 5030.2A at the redshift of RZ 2109). 
The model is in good agreement with the STIS continuum level, 
but there is a clear excess of [O ni] emission. The arrows indicate 
a 10 pixel and 20 pixel wide region, centered on rest frame 5007A. 

& Proffitt 2011). Applying this correction produced the 
final flux calibrated spectrum. 

Figure 1 shows the [O ill] A5007 region of the extracted 
STIS spectrum of RZ 2109. In this figure, we have also 
plotted a previously obtained Keck/LRIS spectrum of 
this cluster (Zepf et al. 2008) and a simple stellar popu- 
lation model from Maraston (2005) of a 12 Gyr cluster, 
with V=21.0 mag and [Fe/H]=-1.3 (similar to the param- 
eters measured for RZ 2109; Steele et al. 2011). The flux 
of the continuum level is in very good agreement with the 
expectation from the cluster model, suggesting that the 
majority of the cluster luminosity is included in the STIS 
slit. The continuum of the Keck spectrum is slightly 
lower at shorter wavelengths. However, this spectrum 
is not fully flux calibrated and was only scaled to have 
the correct flux at 5007A, based on its V-band magni- 
tude (Zepf et al. 2008). It is therefore in good agreement 
with the STIS observations. Gompared with the spec- 
trum expected for a stellar population with RZ 2109's 
properties (black line), it can be seen in both the old 
(Keck, blue line) and new (STIS, red line) spectra, that 
there is strong and broad excess emission around 5007A 
(5030.2A at the velocity of RZ 2109). A discussion of 
the line strengths and any potential temporal variabil- 
ity in these spectra will be presented in Steele et al. (in 
preparation). 

3. STELLAR PROFILE OF RZ 2109 

The size, density, and mass of a GG can significantly in- 
fiuence its population of accreting compact objects. For 
example, it is now well established that the formation 
of LMXBs is orders of magnitude more efficient in GGs 
than the field of the galaxy (e.g. Glark 1975; Angelini 
et al. 2001; Sarazin et al. 2003; Kim et al. 2006; Kundu 
et al. 2007). This is likely due to dynamical interactions 
enhancing the formation of LMXBs. Therefore, clusters 
with higher stellar interaction rates are expected to host 
more LMXBs, as observed (Verbunt & Hut 1987; Jordan 
et al. 2007; Peacock et al. 2009). The presence of an inter- 
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mediate mass BH (IMBH), or even stellar mass BHs, in 
a GC can also significantly affect the cluster's structure, 
potentially producing relatively large cores with shallow 
central cusps (e.g. Hurley 2007; Trenti et al. 2007). Given 
RZ 2109's prominence as one of the few clusters known 
to host a BH, it is interesting to consider the structure 
of the cluster. 

Due to the small angular size of GCs, it is observation- 
ally very challenging to measure the structure of clus- 
ters at the distance of NGC 4472 (which is taken to be 
16.0 Mpc throughout this paper). However, it is possi- 
ble to estimate such parameters by fitting PSF convolved 
King models to high spatial resolution HST images, pro- 
vided the S/N of the observation is high enough (Carl- 
son & Holtzman 2001). Such methods have previously 
been used to estimate the structure of GCs around other 
Virgo cluster galaxies (e.g. Jordan et al. 2004; Waters 
2007; WiUiams et al. 2007; Madrid et al. 2009). 

The target acquisition images that were taken prior to 
each STIS spectrum provide some of the highest S/N 
and spatial resolution images available for this clus- 
ter. In total eight 240 s exposures were obtained on 
the STIS CCD through the F28X50LP fiher (covering 
5500A < A < lOOOOA). The STIS pixels are slightly un- 
dersampled. We therefore drizzled the images together 
onto a grid with pixels 1/10 the native pixel size. Because 
the orientation and location of the source on the detec- 
tor varied between each exposure, this process reduces 
pixellation effects in the final images. Before combina- 
tion, all images were centered by using the iraf/stsdas 
task ELLIPSE to determine the center of the cluster in 
each exposure. The resulting image has a total expo- 
sure time of 1920 s. The cluster profile, obtained from 
this combined image using ellipse, is plotted in Figure 
2 (black plusses). 

We compared the cluster profile to HST/WFPC2 ob- 
servations of the cluster through the F555W filter (pro- 
posal ID 11209, PI Zcpf). Waters et al. (in preparation) 
have previously fit PSF convolved King models to this 
WFPC2 image of the cluster using the Waters (2007) 
code SUPERKING. The resulting structural parameters of 
RZ 2109 are listed in Table 1. This King model, con- 
volved with the PSF for these STIS observations (pro- 
duced by the TinyTim PSF simulator), is plotted in Fig- 
ure 2. It can be seen that the model from the analysis of 
the WFPC2 images provides an excellent representation 
of the STIS cluster profile. Although a full discussion of 
the structure of this cluster is beyond the scope of this pa- 
per, the agreement of the WFPC2 model fit to the STIS 
data is reassuring. This is because the two datasets have 
different systematic uncertainties due effects such as PSF 
shape and pixellation, which typically dominate the un- 
certainties in the fits. Thus the agreement between the 
WFPC2 model and the STIS data provides some assur- 
ance that the parameters are reasonably well determined. 
In this context, we note that RZ 2109 appears to have a 
high concentration and a large half light and tidal radii, 
as may be expected given its large distance from the cen- 
ter of the host galaxy (>30 kpc). 

4. SPATIAL DISTRIBUTION OF THE [O III] NEBULA 

At the distance of NGC 4472 (16 Mpc), RZ 2109 is spa- 
tially resolved by the HST STIS (with 1 pixel - 3.9 pc). 
The resulting long slit spectrum therefore allows us to 
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Fig. 2.— Radial profile of RZ 2109, as measured from the STIS 
target acquisition images (black pluses). The dotted red line in- 
dicates the point spread function (PSF) of these observations and 
the solid blue line shows the previously determined King model fit 
to RZ 2109, convolved with the PSF. This King model fit provides 
an excellent match the new STIS images. 



TABLE 1 
Structural parameters of RZ 2109 



My' Wj; c'" ri," (pc) rj" (pc) rj" (pc) 
-9.76 8.4 1.9 0.92 6.49 82.2 

Derived from King model fits to WFPC2 observations of 
RZ 2109 (Waters et al., in preparation). Parameters are the 
'absolute V-band magnitude, "central potential, "'concentration 
c = log{rt /r^), '^core, half light and tidal radii of the cluster. 

study the cluster's spectrum as a function of radius 
from its center. Of particular interest is the size of 
the [O III] emitting nebula. From Figure 1, we esti- 
mate that the emission from the [O III] doublet influences 
the cluster spectrum over a broad wavelength range of 
4940 < A[oiii] < 5070A. At other wavelengths, the spec- 
trum is consistent with that expected of a GC with the 
luminosity and metallicity of RZ 2109 (Maraston 2005; 
Steele et al. 2011). In Figure 3, we plot the width of 
the spectrum (averaged over 50 pixels = 136 A) in four 
regions centered at 4756, 4893, 5168 and 5305A. These 
regions are at wavelengths above and below the observed 
nebular emission and therefore represent emission from 
only the stellar component of the cluster. We note that, 
while the PSF is known to increase with increasing wave- 
length, the PSF and cluster profile do not vary signifi- 
cantly over this relatively narrow wavelength range. The 
dashed red line in this figure is the profile of a stellar 
source, the T/S'T standard star BD-I-75D325. This spec- 
trum was obtained under a similar setup to RZ 2109's 
and is representative of the PSF of our observations. 
The solid blue line is the King model profile for this clus- 
ter, produced from fits to WFPC2/STIS photometry and 
convolved with the stellar PSF. The cluster is clearly re- 
solved by these observations and the profile of the spec- 
trum is consistent with the previously determined King 
model representation of the cluster. 

To investigate the nebula in RZ 2109, we focus on a 
20 pixel (55A) region centered on the main 5007A peak. 
This region includes nearly all of the 5007A fiux (see Fig- 
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Fig. 3. — The spatial profile of RZ 2109's spectrum across 50 pixel 
wide regions centered at 4756, 4893, 5168 and 5305A (blue crosses). 
These regions are outside of the wavelength range influenced by the 
nebular emission and represent the cluster's continuum emission 
from its stellar component. The blue line represents a fit to the 
average profile across all four of these continuum regions. The 
black points show the average profile of a 20 pixel region centered 
on the [O III] emission line. The nebular emission is centered on 
the cluster and very bright, increasing the cluster flux over these 
wavelengths by around 40%. We have also plotted the estimated 
PSF for these observations (scaled to match the peak flux, dashed 
red line) and PSF convolved King model of RZ 2109 (solid blue 
line). It can be seen that the cluster is clearly extended and is 
well represented by the previously determined King model of the 
cluster. In figure 4, we compare the profile at 5030A to this PSF, 
King model and other profiles. 

ure 1). We also studied 27 A and 14 A regions, which con- 
tain less of the broad [O ill] emission, but have a higher 
ratio of line to cluster flux. The profiles for the differ- 
ent width regions were found to be consistent. In Figure 
4, we plot the profile of this nebular emission. This is 
the average flux over the 20 pixel [O ill] region minus 
the cluster stellar component, taken to be the average of 
the continuum 100 pixels either side of the line region. 
The cluster component at these wavelengths is expected 
to be relatively flat and have no emission or absorption 
lines with significant flux (e.g., Maraston 2005). There- 
fore, this subtracted proflle should accurately represent 
the nebular emission. In this figure, we compare the pro- 
file with several model profiles. All of these models have 
been scaled to have the same integrated flux as the data. 
The reduced values from fitting a selection of these 
flux profiles to the data are listed in Tabel 2. 

It can be seen that the nebula is best represented as 
an extended source. The top left panel of Figure 4 shows 
that the point source representation (dashed red line) has 
a narrower peak than the data and provides a poor fit 
(with a reduced = 4.8, rejecting the fit with a prob- 
ability 99.99%). The solid blue line in this panel shows 
a PSF convolved King model, which represents the clus- 
ter's stellar component. The nebular emission is quite 
well represented by this King model, suggesting that the 
nebula has a similar scale to the cluster's stellar compo- 
nent (i.e., a half light radius of 6.5 pc). However, we do 
not know the functional form of the [O ill] spatial distri- 
bution a priori, therefore we consider several other forms 



TABLE 2 

Reduced values for a selection of flux profiles 



Model ^ 


Parameters 


^ 2 /. . ii 


Stellar psf 




4.81 


Gaussian 


hwhm=2.U pc 


2.20 




hwhm=2.5 pc 


1.55 




hwhm=3.0 pc 


1.13 




hwhm=3.5 pc 


0.94 




hwhm=4.0 pc 


0.96 




hwhm=4.5 pc 


1.17 




hwhm=5.0 pc 


1.52 




hwhm=5.5 pc 


1.99 




hwhm=6.0 pc 


2.55 




hwhm=6.5 pc 


3.17 


r QC L/T 


U-i-cr-c D.O pc 


z.zo 




0.1<r<20.0 pc 


1.02 




0.1<r<80.0 pc 


1.23 




0.5<r< 6.5 pc 


1.82 




0.5<r<20.0 pc 


1.33 




0.5<r<80.0 pc 


2.69 




1.0<r< 6.5 pc 


1.26 




1.0<r<20.0 pc 


1.92 




1.0<r<80.0 pc 


4.11 


King model 


as in Table 1 


1.58 



'AH models for the emission profile have been convolved with 
the PSF for these observations. "It should be noted that only 
the background level is fit during these comparisons (the flux is 
flxed to that observed and the parameters defining the shapes 
of the profiles are fixed at the values quoted). We therefore 
take the number of degrees of freedom for each model, u = 6. 
"'Representing the stellar profile of RZ 2109. 

for the radial profile. In the top right panel of Figure 
4, we compare the nebular emission to PSF convolved 
Gaussian functions of different widths. As can be seen 
from Table 2, the emission is best represented by Gaus- 
sians that have a half width half maximum (HWHM) 
of 4±1.5 pc. We also consider, in the bottom panels 
of Figure 4, the profile produced if the nebula emission 
followed an inverse square law from the center of the clus- 
ter. As an example, we consider three groups of models 
which have density profiles which truncate at 6.5, 20.0 
and 80 pc. For each of these models, we start the emis- 
sion at three different inner radii (0.1, 0.5, 1.0 pc). The 
models truncated at large radii still have significant flux 
at larger distances. While our data show no evidence 
for such emission, the errors on the measured fluxes do 
not allow us to discard such large scale faint emission. It 
can be seen that inverse square law models which have 
small inner radii and large outer radii provide a reason- 
able fit the observations. Similarly, the observations can 
be represented by smaller outer radii if the emission is 
suppressed, by some mechanism, to relatively large inner 
radii. Due to this degeneracy, we are unable to determine 
the exact morphology of the cluster -with several of the 
models providing a good fit. However, all of the models 
which provide a good fit to the data have half light radii 
'--^3-7 pc. 

The actual nebula may be more complex than the mod- 
els considered (e.g., Steele et al. 2011). If it were pro- 
duced by winds from a central BH-LMXB (as has been 
proposed, e.g., Zepf et al. 2007), then these winds may 
not produce a spherically symmetric gas distribution. 
Simflarly LMXBs, particularly BH LMXBs (King et al. 
1996), can have highly variable accretion rates. The gas 
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Fig. 4. — Profile of the [O III] emission line region (black points). The difi'erent panels compare the observed profile with those expected 
from: a stellar source (dashed red line, top left); the PSF convolved King model representing RZ 2109's stellar component (blue line, top 
left); PSF convolved Gaussians (top right); and profiles with varying inner and outer boundaries (bottom panels). All of the models 

are scaled to have the same integrated flux as the data. The reduced values from fitting these profiles to the data are listed in Table 2. 



expelled during such accretion is likely proportional to 
the accretion rate, potentially producing a nebula with 
highly variable density as a function of radius. The data 
do not allow us to investigate such details, so we do not 
consider more complicated models. Instead, we focus on 
the basic constraints provided by these observations - 
that the nebula is clearly extended on parsec scales and 
best represented by having half of its emission on scales 
of -5±2 pc. 

5. THE [O III] NEBULA AND NATURE OF RZ 2109'S BH 

Since the discovery of the ultraluminous X-ray source 
and broad [O ill] emission in RZ 2109, several processes 
have been suggested as the possible source. It is thought 
that the strong and variable X-ray emission requires an 
accreting BH origin. However, the nature of this BH 



hole remains uncertain, with both stellar mass BHs and 
IMBHs proposed. While all current models propose that 
the observed nebula is ionized by the central accreting 
source, several models of the distribution of the nebular 
gas have been proposed. The spatial scale of the [O ill] 
emission, as measured from these observations, can help 
to distinguish between these models. 

5.1. Ionization of a wind from a stellar mass BH 

The [O III] emission observed in RZ 2109 may result 
from the photoionization of a wind driven from the cen- 
tral accreting object (Zepf et al. 2008). BHs, accret- 
ing at rates near the Eddington limit, are thought to 
drive strong winds from the systems (e.g., Proga 2007; 
Zepf et al. 2008, and references therein). Such a wind 
would then be photoionized by the extremely bright ra- 
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diation that is observed from the central source. Zepf 
et al. (2008) demonstrated that this model is capable of 
producing both the luminosity and width of the observed 
[O III] emission. However, driving this strong wind from 
the system requires that the compact object be accret- 
ing at a rate similar to its Eddington limit. Therefore, 
at the X-ray luminosities observed, this model predicts 
accretion on to a stellar mass BH. 

The velocity of the outflowing wind is measured to be 
~ 10^ kms~^, which is far greater than the escape ve- 
locity of a typical cluster of this size, a few xlO kms~^. 
Therefore, the size of the interstellar gas is simply given 
by this wind velocity times the time that the X-ray bi- 
nary has been accreting and driving the wind. Unfor- 
tunately, this timescale is poorly constrained. We can 
estimate a minimal timescale, based on X-ray detections 
of this source during outburst. The source was detected 
by ROSAT in June/ July 1996 , XMM-Newton in June 
2002 and January 2004, marginally by Chandra in Febru- 
ary 2010 (a complete list of X-ray detections of RZ 2109 
prior to 2010 is presented in Maccarone et al. 2010) 
and was still detected in the last observation in 2011 
(Tana Joseph, private communication). We therefore as- 
sume that the source has been in outburst throughtout 
this period, with its obsevered luminosity peaking at 
^4xl0''^erg/s and reducing to '--^IxlO'^^erg/s. Assum- 
ing steady accretion over 14 yrs would suggest that the 
emitted gas has a minimum size of ~0.05 pc. Setting an 
upper limit is much harder because the evolution of X-ray 
binaries is highly variable and uncertain. Some sources 
can be persistent for long timescales, while others are 
highly variable and switch between periods of quiescence 
and outburst. XMMU122939.7-h075333's history prior 
to the i?05'^T' observation is therefore unknown. Given 
the broad range of possible accretion (and hence wind) 
timescales, this model could have resulted in gaseous re- 
gions at a range of radii, including those observed here. 

5.2. X-ray ionized nova ejecta 

Another potential source of the nebular gas is the 
ejecta from a nova in the cluster (Ripamonti & Mapelli 
2012). In this model, the accreting object photoionizes 
the ejecta from an unrelated nova in the same cluster. 
Because the ionizing source and nebula gas are unrelated 
in this model, the accreting compact object may have 
any mass, so long as it produces sufficient X-ray lumi- 
nosity. Ripamonti & Mapelli (2012) demonstrated that 
this model is capable of reproducing both the [O in] dou- 
blet and the lack of other emission lines in the observed 
optical spectrum of RZ 2109. However, their model re- 
quires that the nova ejecta shell is <0.1 pc from the X-ray 
source. This is significantly smaller than the observed 
nebula, which is found to be on parsec scales. 

5.3. IMBH accretion disk origin 

Previous studies have considered the possibility that 
the large velocity width of the [O iii] A5007 emission is 
due to gravitational motions around a BH. Zepf et al. 
(2008) foimd that such a model fails because, at the 
radius where such velocities are available, there is not 
enough volume to prodiice the observed [O III] luminos- 
ity (given its critical density) unless the mass of the BH 
is a significant fraction of the total GC mass. Porter 



(2010) confirmed this calculation, and showed specifically 
that such a scenario was only possible with a BH mass 
>30,000 Mq. If the width of the [O ill] emission is due 
to Keplerian motions around a 30,000 BH, then the 
emitting region would be around 10~^ pc. Such a region 
is many orders of magnitude beneath the spatial reso- 
lution of our observations and would appear unresolved. 
Given that the [O ill] emission line is resolved, these ob- 
servations rule out such an origin. 

5.4. Tidal disruption of a white dwarf/ horizontal 
branch star by an IMBH 

Another model, invoked to explain ultraluminous X- 
ray emission and optical emission lines, is the tidal dis- 
ruption of a white dwarf by a central IMBH (e.g., Irwin 
et al. 2010). In this model, debris from the tidally dis- 
rupted star is accreted on to the central IMBH. The re- 
sulting outburst can then photoionize the remaining de- 
bris, hence producing both the observed X-ray and neb- 
ular emission. Clausen & Eracleous (2011) noted that, 
while such a model can produce both the X-ray emission 
and the luminosity and width of the [O ill] A5007 emis- 
sion line observed from RZ 2109, the emission should 
peak after a few years. Such a timescale is not long 
enough to explain the observed X-ray emission from the 
cluster, which suggests that the source has been in out- 
burst for over 14 yrs. It is also inconsistent with the neb- 
ular emission observed here, whose scale of > 2 pc and 
velocity of 1500 kms~^ suggest a timescale of > 1300 yrs. 

Clausen et al. (2012) extended this work to consider the 
tidal disruption of a horizontal branch star by an IMBH. 
They demonstrated that this could produce the X-ray, 
[O III] and [N ii] emission observed from another GC BH 
candidate, the X-ray source CXOJ033831. 8352604 in an 
NGC 1399 GC. The disruption of a horizontal branch 
star predicts a longer outburst timescale with r ^200 yrs 
(Clausen et al. 2012) and is therefore capable of produc- 
ing the observed X-ray emission. At the observed veloc- 
ity of the emission line in RZ 2109, a ^1500 kms"-'^ (Zepf 
et al. 2008), the maximum associated nebula scale is sim- 
ply T X CT < 0.3 pc. Such a nebula would be unresolved 
by these observations and is smaller than the best fitting 
half light radii to the observed [O III] emission. There- 
fore, the tidal disruption of a star by an IMBH is unlikely 
to be the source of the emission observed in RZ 2109. 

6. CONCLUSIONS 

The STIS spectroscopy presented in this paper spa- 
tially resolves both the stellar and nebular emis- 
sion from the GC RZ 2109. By considering the 
[O III] AA4959, 5007A emission lines, we have shown that 
the nebula in this cluster is extended, with best fitting 
half light radii in the range 3-7 pc. 

The observed size of this nebula is larger than that 
predicted by several models that have been proposed to 
explain the nebular and ULX emission detected from this 
cluster. It is many orders of magnitude larger than that 
expected from emission from an accretion disk around an 
IMBH (which should be on scales of lO^'' pc). Similarly, 
the tidal disruption of a WD or HB star by an IMBH pre- 
dicts emission on timescales that produce nebulae that 
are an order of magnitude smaller than observed (0.05 pc 
and 0.3 pc, respectively). The nebula is also an order of 
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magnitude larger than that predicted from the ioniza- 
tion of ejecta from a nova in the cluster, which should 
produce nebulae on scales of 0.1 pc. 

The large nebula observed could be produced via a 
wind driven from a BH accreting at a rate close to the 
Eddington limit. The timescales for such accretion are 
unknown. However, it is possible that such a system 
has been accreting for long enough (either persistently 
or with transient behavior) to drive the nebula gas to 
the few pc scales observed. 
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